Random Walks and Anderson Localisation in a Three-Dimensional Class C Network Model 
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We study the disorder-induced localisation transition in a three-dimensional network model that belongs to 
symmetry class C. The model represents quasiparticle dynamics in a gapless spin-singlet superconductor without 
time-reversal invariance. It is a special feature of network models with this symmetry that the conductance and 
density of states can be expressed as averages in a classical system of dense, interacting random walks. Using 
this mapping, we present a more precise numerical study of critical behaviour at an Anderson transition than 
has been possible previously in any context. 
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Anderson transitions between diffusive and localised 
phases of quantum particles in disordered systems constitute 
an important category of critical phenomena. Depending on 
dimensionality and the symmetries of the Hamiltonian, vari- 
ous universality classes are possible for scaling behaviour |TJ. 
In most cases there is no known connection between these 
Anderson universality classes and those for phase transitions 
in classical systems, but in a special instance, known as class 
C, properties of suitably chosen quantum lattice models can 
be expressed in terms of observables for a classical model de- 
fined on the same lattice. This mapping was originally discov- 
ered in the context of the spin quantum Hall effect, where it re- 
lates a delocalisation transition in two dimensions to classical 
percolation, also in two dimensions, for which many relevant 
aspects of critical behaviour are known exactly [2j. In this 
paper we apply the mapping to a three-dimensional system 
for which the classical counterpart is a model of interacting 
random walks (3). The classical model is of interest both as 
a representation of the quantum problem and in its own right. 
While no exact results for its behaviour are available, the map- 
ping makes possible simulations of the Anderson transition 
with unprecedented precision. 

The symmetry class for localisation that we are concerned 
with is one of seven recognised about a decade ago as be- 
ing additional to the three Wigner-Dyson classes that were 
originally identified in the context of random matrix theory 
|4|. Systems in these additional symmetry classes are distin- 
guished from ones in the Wigner-Dyson classes by having a 
special point in the energy spectrum and energy levels that ap- 
pear in pairs either side of this point. In particular, models in 
class C arise from the Bogoliubov de-Gennes Hamiltonian for 
quasiparticles in a gapless, disordered spin-singlet supercon- 
ductor with broken time-reversal symmetry for orbital motion 
but negligible Zeeman splitting. Here the special energy is 
the chemical potential (which we set to zero) and pairs of lev- 
els are related by particle-hole symmetry. Quasiparticle states 
in a three-dimensional system of this type have an Anderson 
transition as a function of disorder strength, which must be 
probed via spin or thermal transport since charge transport is 



short-circuited by the condensate. Particle-hole symmetry has 
profound consequences for the influence of disorder on quasi- 
particle eigenstates and for the Anderson transition, which 
have been explored previously using random matrix theory 
1 4 1, the non-linear sigma model 0161, and via calculations for 
one-dimensional systems [7]. Most importantly, the density of 
states at the chemical potential plays the role of an order pa- 
rameter for the transition, being finite in the metal and zero 
in the insulator; by contrast, the density of states shows no 
critical behaviour at transitions in the Wigner-Dyson classes. 

The quantum to classical mapping provides a framework 
within which quasiparticle properties can be studied in great 
detail starting from a simplified description of a disordered 
superconductor. It is based on a formulation of the localisa- 
tion problem as a network model |8| in which quasiparticles 
propagate along the directed links of a lattice and scatter be- 
tween links at nodes. Disorder enters the model in the form of 
quenched random phase shifts associated with propagation on 
links. For versions of the network model belonging to class C, 
first introduced in Ref. [9|, the disorder-averaged quasiparti- 
cle density of states and (spin) conductance can be expressed 
as averages over configurations of interacting classical ran- 
dom walks on the same directed lattice J2][3]. This relation 
between quantum properties and classical walks holds on any 
graph in which all nodes have exactly two incoming and two 
outgoing links. A single parameter p controls behaviour at 
nodes: incoming and outgoing links are arranged in pairs, 
and a particle passing through the node follows the pairing 
with probability p, or switches with probability 1 — p. More 
precisely, in the quantum problem the probabilities calculated 
by squaring amplitudes in the scattering matrix for the node 
take the values p and 1 — p, while in the classical problem the 
connection at each node between ingoing and outgoing links 
is a quenched random variable having two possible arrange- 
ments with these probabilities. Given a directed graph with 
the required coordination, any choice of classical connections 
at the nodes separates paths on the graph into a set of distinct, 
closed, mutually avoiding walks. Average properties of these 
walks are calculated from a sum over all node configurations, 



2 



weighted according to their probabilities. 

To use the mapping for a three-dimensional model one must 
select a four-fold coordinated lattice, assign directions to the 
links so that two are incoming and two are outgoing at every 
node, and pair incoming with outgoing links. We pick the dia- 
mond lattice and choose link directions and pairings at nodes 
so that the system has localised states at p = and extended 
states at p — 1. The resulting unit cell contains 24 sites and we 
defer a full descrition until later in this paper. We study sam- 
ples consisting of L x L x L unit cells with two alternative 
sets of boundary conditions: using the classical mapping, we 
calculate conductance between two opposite, open faces (with 
periodic boundary conditions in the other directions), and we 
calculate the density of states in closed samples with periodic 
boundary conditions in all three directions. 

The results of the mapping are as follows |3|. First, the 
disorder-averaged spin conductance G(p, L) of the quantum 
system is given (in units of H/An) by the average of the num- 
ber ATl of classical paths from a specified open face to the 
other. Second, the average density of states p(e) for eigen- 
phases of a unitary time evolution operator, which occupy the 
range — 7r < e < ir and play the role of energy levels, is given 
in terms of the probability P(s,p) in the classical problem for 
a given link to belong to a closed walk of length s steps, by 
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^ P(s, p) cos(2se) 



s>0 



(1) 



In the following we present and analyse data for G(p, L) and 
P(s,p), and deduce from the latter the critical properties of 
p(e, p). In addition, we study the mean square end-to-end dis- 
tance (R(s) 2 ) of walks as a function of number of steps s. 
This last quantity has no direct significance for the original 
quantum problem but provides a useful characterisation of the 
classical walks. 

We can study much larger systems using the classical de- 
scription than is possible for a quantum Hamiltonian with dis- 
order: we use sizes in the range 40 < L < 440, and calculate 
averages over 10 7 realisations for the smallest samples and 
4 x 10 3 for the largest. The biggest systems hence contain 
over 2 • 10 9 sites in total, and over 10 6 sites in cross-section. 
By contrast, direct calculations of eigenfunctions using sparse 
matrix techniques ifTOl [HI [T2l are limited to systems 10 2 or 
10 3 times smaller, while high precision studies of Anderson 
transitions lfl3l using transfer matrix methods lfl~4l are re- 
stricted to systems with a cross-section of around 10 3 sites. 
Our system sizes enable us to determine the values of critical 
exponents with about an order of magnitude higher precision 
that has been possible numerically for other Anderson transi- 
tions. 

We discuss first our data for conductance G(p,L). Our 
measurements span the range 0.25 < p < 0.28 and yield 
values for G(p, L) in the metal with statistical uncertainties of 
less than 0.1%, except for the largest sample size (L = 440) 
where they are 0.4%. The model is designed to have a critical 
point at p = p c separating a metallic phase for p > p c from 
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FIG. 1: (Color online) Conductance as a function of (p — p^L 1 ^, 
illustrating scaling collapse. Upper inset: same data on a logarithmic 
conductance scale. Lower inset: conductance as a function of p for 
several values of L. Lines are scaling fit described in main text. 



an insulating one for p < p c . From Ohm's law the metallic 
phase is characterised by the behaviour G(p, L) ~ <r(p)L for 
sufficiently large L, where a(p) is the conductivity, while in 
the insulator G(p, L) — ► at large L. Close to the critical 
point we expect one parameter scaling in the form 



G(p,L) = f(L/£(p)), 



(2) 



where £(p) is the localisation length in the insulator and the 
correlation length in the metal. Hence curves of G(p, L) as a 
function of p for different L intersect at p c , as illustrated in the 
lower inset to Fig. |T| In addition, assuming that diverges 
at the transition as £(p) oc \p — p c \~ h ' , data for G(p,L) in 
the critical regime should collapse onto a single curve when 
plotted as a function of x = (p — p c )L x l v . This collapse 
is shown in Fig. [T] using a linear scale for G(p,L), and in 
the upper left inset using a logarithmic scale to display more 
clearly behaviour in the insulator. 

Our estimates for the values of v and p c obtained from such 
a scaling analysis are 



^ = 0.9985 ±0.0015 



(3) 



and p c = 0.260116 ± 0.000002. The analysis in detail is as 
follows. Extending Eq. Q to include corrections to scaling, 
we expect 



G(p,L) = /(L 1 /^,i^V) 



(4) 



where <fi is the relevant scaling variable and t/j is the scaling 
variable for the leading correction, which is irrelevant pro- 
vided y m < 0. Since G(L,p) varies smoothly withp at fixed 
L, the scaling variables have Taylor expansions 



<t>= (p-Pc) + A{p- Pc f 



(5) 



and ip = 1 + B(p — p c ) . . ., in which leading coefficients have 
been fixed via the definition of f(x, y). As a first step, we omit 
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all corrections to scaling, setting A and higher coefficients to 
zero and considering f(x, 0). We construct this scaling func- 
tion using cubic B-splines. We judge all fits by comparing the 
value of x 2 with the number of degrees of freedom. For the 
splines we use 18 internal knot points at a spacing of 0.1 for 
—0.5 < x < 1.0, where curvature of f(x, 0) is highest, and a 
spacing of 1 for 1 < x < 5, where f(x, 0) is almost linear in 
x. A coarser mesh of knot points does not allow an adequate 
fit to the data, while a finer one does not significantly improve 
the fit. Next we allow for an irrelevant scaling variable, by 
Taylor expanding f(x, y) in y: we find that the form 

f(x,y)=f(x,0)[l + Cy] (6) 

with y 1TI = — 1 and C constant provides an adequate descrip- 
tion of the data. Finally, we allow for a non-linear scaling 
variable, as in Eq. |5]), finding from the effect on x 2 that non- 
zero A is justified but non-zero B or further terms are not. 

Taking data with p < p max = 0.273 (excluded values are 
far from the critical point) and x > —0.5) (excluded values 
are deep in the insulator, where errors in G(L,p) are large), 
we have 262 points to fit and 26 parameters (the values of v, 
p c , A, C, 18 internal spline points and four boundary ones). 
Our fit has \ 2 — 236.4 and is insensitive to the value of p max 
in the range 0.265 < p max < 0.28. Since the fit yields a 
value for \ 2 close to the number of degrees of freedom, 262 — 
26 = 236, we believe that errors in our values for v and p c 
are due mainly to statistical errors in the data. We place 95% 
confidence limits on these values in two independent ways. In 
one we find the variation that increases \ 2 by 4. In the other 
we use a Monte Carlo technique to generate synthetic data 
|[T5l . The two methods give the same results. 

Our value for v is in striking, although presumably acciden- 
tal, agreement with the leading order result from an e expan- 
sion in 2 + e dimensions [5|, v — 1/e, evaluated at e = 1. It 
is also consistent with the value v « 0.9 found in an earlier 
transfer matrix calculation |16| for a three-dimensional class 
C network model, which did not use the mapping to classical 
random walks and so was restricted to much smaller system 
sizes. The closeness of our value to v = 1 is tantalising but we 
are not aware of any theoretical reason to attach significance 
to this. Indeed, the upper critical dimension for the localisa- 
tion transition in class C is believed to be four |[T7ll . and so 
non-trivial exponent values are expected in three dimensions. 

We next present our data for the integrated return probabil- 
ity N(s,p), which is related to P(s,p) via 

N(s,p)=J2 p (t>P) ■ ( 7 ) 

t>s 

From this we deduce the critical behaviour of the classical 
walks and of the quantum density of states. These data are 
obtained in closed samples by generating all walks for a given 
realisation of node configurations and averaging over differ- 
ent realisations, using sample sizes L < 300. In the critical 
regime walks are fractal with dimension df. The length 
sets their characteristic size in the insulator, so that the char- 
acteristic arc length is in the metal £(p) represents 
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FIG. 2: (Color online) \p - p c | (df ~ 3)v N(s, p) asafunction of s\p- 
p c | * " on double logarithmic scales for several values of p either side 
of p c . Inset: N(s, p c ) as a function of s on double logarithmic scales. 

a crossover length beyond which scaling properties are as for 
free random walks. Close to the critical point we expect the 
scaling form 

N(s,p) = Z(p) d <- 3 h ± (s/Z(p) d <) (8) 

to hold, with scaling functions h+(x) for p > p c and h-(x) 
for p < p c that have distinct behaviour at large x but share 
the same power law form, h(x) ~ x 2 ~ T , at small x. Since 
N(s,p) should be independent of £(p) at small s, we have the 
exponent relation 



Data at p = p c are shown in the inset to Fig. [2] they con- 
firm power law behaviour over six decades in s, yielding the 
exponent value 

t = 2.184 ±0.003 (10) 

and hence d{ = 2.534 ± 0.009. Using our values for df and 
v, we demonstrate in the main panel of Fig.|2]scaling collapse 
following Eq. ([8]) for seven values of p in the insulator (0.250, 
0.253, 0.255, 0.256, 0.257, 0.258, and 0.259) and five val- 
ues in the metal (0.261, 0.262, 0.263, 0.265, and 0.270). We 
note the excellent overlap of data in each phase. Both scaling 
functions h±(x) follow the power-law behaviour of the criti- 
cal point at small x. At large x, h-{x) falls off exponentially, 
while h+(x) tends to a constant: the escape probability. 

These results can be combined with those of Ref. to 
conclude for the quantum density of states that: (i) p(e,p) oc 
e 2 in the insulator; (ii) p(e,p c ) oc |e| T_2 at the critical point; 
while in the metal (iii) p(0,p) oc (p — p c )( 3 ~ d f*> l/ and (iv) 
[p(e,p) — p(0,p)] oc | er | 1 / 2 . For comparison, the e expansion 
gives at leading order: p(e,p c ) oc |e| e / 2 and p(0,p) oc 
(p-Pc)- 

To examine further the properties of classical walks at the 
critical point and in the metal we have calculated (R(s) 2 ) for 
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FIG. 3: (Color online) {R 2 (s)) / £(p) 2 as a function of s/£(p) d{ in 
the metal (p = 0.265, 0.27, 0.3 and 0.5), on double logarithmic 
scales. The continuous and dashed lines have gradients of 1 and 
2/df, respectively. Inset: (R(s) 2 ) vs s at the critical point. 



Hexagon 


Initial Sites 


Hexagon 


Initial Sites 


a, —8, f3 


(2,0,2) 


(2,-2,4) 


a,-7,<5 


(0,0,0) 


(0,-2,2) 


5,-~f,P 


(4,2,2) 


(4,4,4) 


«,-/9,7 


(4,2,2) 


(4,4,4) 



TABLE I: hexagons and their initial sites in the unit cell. 



very long trajectories, which can be generated individually in 
samples of unbounded size. We class a trajectory as long if 
it has not closed on itself after 5 x 10 7 steps, and we analyse 
behaviour for the first 10 7 steps. At p = p c we expect to see 
with this method the fractal behaviour of critical walks, and in 
the inset to Fig. [5] we show that (R(s) 2 ) varies as a power of 
s. From this we extract 2/df = 0.784 ± 0.004, which is con- 
sistent with the value determined from t using Eq. |9]). In the 
extended phase we expect scaling collapse of (R(s) 2 ) / £(p) 2 
plotted as a function of s/£(p) dt for different p. This is illus- 
trated in the main panel of Fig. fusing four values of p. Two 
power law regimes are apparent: (R(s) 2 ) ~ s 2 / dt for small 
s/£,(p) dt , as at the critical point; and (R(s) 2 ) ~ s for large 
s/£,(p) dl , as for free random walks. 

We close with a full description of our model, deferred 
above. Consider a diamond lattice. We take as Cartesian coor- 
dinates for the two sites in the primitive unit cell: r = (0,0,0) 
on sublattice A, and r = (1,1,1) on sublattice B. At 
an A-site the position vectors to neighbouring £>-sites are: 
a = (1,1,1), = (-1,-1,1), 7 = (-1,1,-1) and 
S = (1,-1,-1). For p = all classical paths in the model 
are closed hexagons. They can be specified by giving an initial 
site and the first three steps. The unit cell of the model consists 
of such eight hexagons, as listed in Table 1. The model does 
not retain the full point symmetry of the diamond lattice, be- 
ing instead tetragonal. It is, however, not very anisotropic: the 
ratio of conductivities in the two distinct directions is about 
1 . 1 and depends very little on p. 

In conclusion, we have used the mapping to classical walks 



as a means to investigate the three-dimensional class C lo- 
calisation transition with high precision. Possible future ex- 
tensions include the study of two and three particle Green 
functions ||T8l and of lattices with higher coordination num- 
ber |[T9l . Most importantly, an analytical understanding of 
the classical problem, perhaps based on ideas from polymer 
physics |20|, would be very desirable. 
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